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We discuss the propagation of gravitational waves in a recently discussed class of theories 
containing N > 2 metric tensors and a corresponding number of standard model copies. Us- 
ing the formalism of gauge-invariant linear perturbation theory we show that all gravitational 
waves propagate at the speed of light. We then employ the Newman-Penrose formalism to 
show that two to six polarizations of gravitational waves may exist, depending on the pa- 
rameters entering the field equations. This corresponds to E(2) representations N2, N3, III5 
, and 116 . We finally apply our general discussion to a recently presented concrete multimetric 

gravity model and show that it is of class N2, i.e., it allows only two tensor polarizations, as 
it is the case for general relativity. Our results provide the theoretical background for tests 
£SJ ' of multimetric gravity theories using the upcoming gravitational wave experiments. 
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. I. MOTIVATION 

■ This article continues a series of articles [l|-y] discussing gravity theories with N > 2 metric 

> ■ r j 

. tensors g and a corresponding number of standard model copies \P , where each copy of the 

' standard model couples only to its own metric tensor and the interaction between the different 

(N ' 

are particularly interested in theories that exhibit repulsive gravitational forces between different 
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standard model copies is mediated solely by an interaction between the different metrics. We 



standard model copies in the Newtonian limit which are of equal strength compared to the attractive 
gravitational force within each standard model copy. While this is not possible [l| for N = 2, it 



may serve as a potential explanation for the observed late-time acceleration of the universe 
for N > 3 and is consistent with high-precision solar system experiments at the post-Newtonian 
level The theories we consider satisfy the following assumptions: 

(i) The action is of the form 

N 

S = S G [g 1 ,...,g N ] + J2SM[9 1 ,* 1 }, (1) 
1=1 
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where Sg is the gravitational part of the action and Sm is the standard model action. 
This assumption guarantees that each standard model copy Vl^ couples only to its own metric 
tensor g 1 , so that the standard model fields satisfy the same field equations as they would 
in a single-metric theory. It further excludes any non-gravitational interaction between the 
different standard model copies, which implies that the different standard model copies appear 
mutually dark. 

(ii) The gravitational field equations are obtained by variation with respect to the metrics g\ h . . . g~l, 
and so are a set of symmetric two-tensor equations of the form K_ ab = 8nGNT ab . 
This employs the well-known principle of stationary action. It follows that the number of 
field equations equals the number of field components. 

(Hi) The geometry tensor K_ ab contains at most second derivatives of the metrics, which can be 
achieved by a suitable choice of the gravitational action. 

This assumption is one of mathematical simplicity and guarantees a reasonable amount of 
technical control over the partial differential field equations. It will be used to restrict the 
possible terms in the linearized field equations of our theory. 

(iv) The field equations are symmetric with respect to arbitrary permutations of the sectors {g , \P ). 
This is another assumption made for simplicity; it employs the Copernican principle in the 
sense that the same laws of nature should hold within each sector. It also follows that the 
interaction between the different sectors will satisfy Newton's principle that action equals 
reaction for the gravitational forces in the Newtonian limit. 

(v) The vacuum solution is given by a set of flat metrics g l ah = rj ab . 

Cosmological constants are excluded because we are interested in multimetric gravity theories 
in which the accelerating universe is modelled by a repulsive interaction between different 
standard model copies, as we have shown in [2j. Further assuming a simultaneous maximal 
set of Killing symmetries for all metrics yields the stated vacuum solution. 

The aim of this article is to examine the propagation of gravitational waves in the weak field limit 
of multimetric gravity theories satisfying the aforementioned assumptions. Note that in addition to 
the recently discussed theories with N > 3 this also includes bimetric gravity theories such as Q-6]; 
see |?| for an overview of bimetric theories and a discussion of their weak field limits. In particular, 
we aim to calculate two properties of gravitational waves which are expected to be accessible by 
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the upcoming detector experiments. One of these properties is the polarization of gravitational 
waves. Metric gravity theories can be classified by the presence of up to six polarizations in terms 
of representations of the little group E(2), which may be distinguished by measuring the electric 
components Roaop of the Riemann tensor [§[ Q]. Although gravitational waves have not been 
observed yet, the sensitivity of present and future experiments is continuously being improved, and 
it is expected that a sensitivity sufficient for the detection of gravitational waves will be reached 
within the next years, hence providing a valuable instrument for testing gravity theories [lol. 

The other property is the propagation velocity v g , which equals the speed of light in general 
relativity, but may differ significantly in multimetric gravity theories containing massive gravi- 
tons 12j, [l3|]. Theories of this type may be tested by comparing the arrival times of gravitational 
Nation and Hgn, tan distant S ap er aovae t a k in g into account the p„ Ssibil it y of a different 



Shapiro delay for both types of radiation 



13 . Il6j ] . Apart from the direct observation of gravitational 



waves, bounds on v g can also be obtained through indirect observations. An upper bound on c — v g 
is placed by the observation of high-energy cosmic rays: massive particles whose velocity exceeds 
v g should be decelerated due to the emission of gravitational bremsstrahlung [17]. Another bound 
on c — v g can be obtained from pulsar timing: if v g < c, the interaction between electromagnetic 
and gravitational radiation should influence the arrival times of radio signals in a gravitational 
wave background Q 5 Q] • 

A complete calculation of the aforementioned effects such as the Shapiro delay of both light and 
gravitational waves in cosmic gravitational fields or the mutual interaction between both types of 
radiation requires a treatment based on the full non-linear field equations of a concrete multimetric 
gravity theory. Since it is the aim of this article to derive the general properties of gravitational 
waves for a large class of multimetric gravity theories, we will not perform this calculation here. 
Instead, we will make use of assumption and consider the propagation of gravitational waves 
in a flat background metric. We then compare their propagation velocity v g to the fundamental 
velocity c, which equals the physical speed of light in this background according to assumption @). 

The outline of this article is as follows. In section [TT] we will determine the propagation velocity 
of gravitational waves in a flat Minkowski background. For this purpose we will apply the gauge 



invariant linear perturbation formalism known from cosmological perturbation analysis [2CH22| 
to the most general linearized vacuum field equations of multimetric gravity in subsection III Ai 
This formalism allows us to separate the physical degrees of freedom from pure gauge quantities. 
In subsection III Bl we will discuss the role of the Bianchi identities in the class of theories we 
consider. We will then calculate the wave-like solutions of the linearized vacuum field equations in 
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subsection III CI and show that all wave-like solutions are null waves. In section Hill we will use the 
Newman-Penrose formalism 23J to determine the allowed polarizations. The findings of sections [III 



and IIIII will then be applied to two concrete example theories in section IIVI We will conclude with 
a discussion in section [V] 



II. PROPAGATION VELOCITY 



In this section we will calculate the propagation velocity v g of gravitational waves within the 
class of multimetric theories satisfying assumptions |§) to (fjuj) stated in the introduction. The 
starting point of our calculation will be a perturbation ansatz around the vacuum solution, which 
is a set of flat Minkowski metrics according to assumption (dj). This ansatz leads us to the most 
general linearized multimetric vacuum field equations satisfying our assumptions. We will employ 
the gauge-invariant formalism detailed in [l] in order to determine the physical degrees of freedom. 
It will turn out that the only wave-like solutions of the gauge-invariant field equations are null 
waves. 

From assumption fz]) it follows further that light rays constituted by the electromagnetic field 
of one standard model copy iff 1 follow the lightlike geodesies of the corresponding metric g^ b . The 
fact that we use the aforementioned perturbation ansatz allows us to conclude that these geodesies 
are the lightlike directions of the Minkowski background, up to higher order perturbations which 
we neglect. Hence, we will conclude that all gravitational waves propagate at the speed of light. 



A. Gauge-invariant formalism 

For the derivation presented in this section it is sufficient to treat gravitational waves as a 
small perturbation of the metrics around a vacuum solution of the field equations. Making use of 
assumption (fgj, we thus use the perturbation ansatz 

9L = Vab + hi b , (2) 

where the components h} ab are small, \h^ b \ >C 1. Under this condition the most general linearized 
vacuum field equations compatible with our assumptions stated in the introduction take 

the form [1] 

o = = e ■ d p d (a h b)p + Q-uh ab + g- d a d b h + m • d p d q h pqVab + k ■ ah Vab , (3) 
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where indices are raised with the flat metric r\ and □ = i] ab d a db. The matrices P, Q, R, M, N 
are constant parameters. Note that these parameter matrices are not completely arbitrary within 
the class of theories we consider, but are further restricted by our assumptions, as we will show 
in the remainder of this section. It will turn out that our assumptions do not uniquely fix the 
parameter matrices. For any concrete multimetric gravity theory, their values depend on the 
choice of the gravitational action Sg introduced in equation (pQ), and can be calculated by an 
explicit linearization of the full nonlinear field equations. We will give the values of the parameter 
matrices for two example theories in section IIVI 

We now apply th e g au g e invariant linear perturbation formalism known from cosmological 



perturbation analysis 



204221] to the linearized field equations ([3]) in order to determine the physical 



l| for full detail. First, we perform a 



degrees of freedom. We only sketch the procedure here; see 
purely algebraic (1 + 3) split of the spacetime coordinates x a = (x°,x a ) into time and space, and 
correspondingly decompose the metric perturbations h ab and the geometry tensors K ab . Second, 
we perform a differential decomposition of the metric perturbations, 

ZLoo = -20, h Qa = d a B + B a , h a p = -2l5 a p + 2A a pE + 4d (a Ep ) +2E a p, (4) 

into four scalars (f>,tp,B_,E, two divergence- free (or transverse) vectors B_ a ,E a and one divergence- 
free, trace-free tensor E^g. Here = d a dg — ^5 a [sA denotes the trace-free second derivative. A 
similar decomposition of the geometry tensor, 



^oo> K 0a = d a W + W a , K ap = \z5 ap + A aB Z + 2d {a Z 0) + , (5) 



shows that its scalar components K_ 00 ,W_, Z_, Z_ depend only on scalar components of the metrics, 
its transverse vector components W_ a ,Z_ a depend only on vector components, and its transverse 
trace-free tensor components Z_ a a depend only on tensor components. In other words, the scalar, 
vector and tensor components of the field equations decouple. In the next step we replace the 
components of the metric perturbations by the potentials 



h= (p + d R- OqE, I 2 = ^ + ^AE, h = B, h = E 



1 

.... — L 

— - 3 

L a = B a - 2d E a , JLu= Ma, La/3 = Ea/3 ■ ( 6 ) 
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Using these quantities, we finally obtain the scalar equations 

^oo = 2(E + g + R + K + ■ d oh~ 2(g + K) ■ ^h- ^{R + K) ■ d oh 

+ 2{K + 3N) • Al2 + 2(E + Q + R + K + K) ■ {-dlh + d^U) (7a) 
+ (P + 2Q + 2^ + 2^) • d A/3 + 2(R - Q) • d$ A/4 - 2(M + • AA/4 , 
W = (£ + 2£) ■ d h - (£ + 6£) • d h - ^(3£ + 2£ + 4£) • dlh 

+ \{P + 2£) • A/3 + (P + 2£) • a (5 2 + A)A , (7b) 
Z = -6(M + £D • ^q/i + 2(g + 3j£) • A/i + 6(Q + 3K) ■ + 2(£ - Q) • ^ A/4 

- 2(£ + 3£ + 3£ + 3M + 9^) • A/2 + 6{M + #) • (d^/3 - d^h) (7c) 

- (£ + 2£ + 6M + • d A/3 + 2(£ + Q + £ + 3M + 3^) • A A/4 , 

Z = 2£ • /1 - 2(£ + 3£) /a - (£ + 2£) ■d h + 2(R-Q) ■ d%h + 2(E + g + £) • A/4, (7d) 

the vector equations 

$U = + 2 £) • ( a o4 - 25 nr a ) + £ • a/ q , (8a) 

la = ~\g ■ 9oLa + (P + 2£) • , (8b) 

and the tensor equations 

l Q/3 = 2£.D/ a/3 . (9) 

The reason for this rewriting becomes apparent when we determine the physical degrees of freedom. 
These are linear combinations of the components of the metrics which are invariant under gauge 
transformations. Since the gravitational fields of our theory are a set of metric tensors, the only 
gauge transformations we allow are diffeomorphisms of the underlying manifold, as it is also the 
case in general relativity. Any such gauge transformation is generated by a vector field £ and 
simultaneously changes all tensor fields F by their Lie derivatives, 5^F = C^F. Considering the 
special case F = g 1 , we find that the components of the metric perturbations transform according 
to 

k h L = datb + d b i a , (10) 

where £ a = rja^ . Note that this gauge freedom is more restrictive than it would be for a set of iV 
independent spin-2 fields, where each field has its own set of gauge transformations [3,0]. Since 
we further wish to keep the formal structure of the perturbation ansatz ([2]), we only consider gauge 
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transformations in which the components £ a are of order 0(h). Applying the decomposition from 
equation (H|) to the vector field £ a , 

& = i* = dd+L, (if) 

we obtain two scalars £ and £ and one divergence- free vector £ a . These can be used to write the 
change of the potentials ([6]) under gauge transformations in the form 

6tl( = 6<:lZ = 0, 8zlZ = do£+S, <5^ 4 7 = |, <^ = 0, S^ = ^ a , 5^ = 0. (12) 

One can now immediately read off the gauge- invariant quantities I\, I2, L a an( l L a s- Further, 
linear combinations of the form C7/3, Cf/J and cjI'J (where summation over I is implied) are 
gauge-invariant if and only if the sum of the coefficients cj vanishes. It thus makes sense to 
consider the linearly related quantities 3 = U_ ■ J, where the matrix U_ is given by 



U 



1.1 



4= if I = 1 or J = 1 

JN 



1 + 4- + 1 if I = J > 1 (13) 
otherwise. 



/AT-AT 

One can easily check that the gauge- invariant degrees of freedom are then given by ^{^i^a^aB 
and 3\, J|, 3'*, where uppercase indices /, J, . . . run from 1 to N, while lowercase indices i,j, . . . run 
from 2 to N. The remaining quantities ^3,^4, 3'^ are pure gauge degrees of freedom and correspond 
to the two scalars and the transverse vector component of the diffeomorphism vector field £. 

The choice of the basis transformation (|13p has another advantage. From assumption it 
follows that the field equations must be invariant under arbitrary permutations of the sectors. For 
the linearized vacuum field equations this implies that the parameter matrices must be invariant 
under simultaneous permutations of both indices, i.e., under transformations of the form 

U ^ o KL a I K <y J L , (14) 

for arbitrary permutation matrices a. It then follows that the entries IJ must be independent of 
the individual values of the indices /, J, and they may only depend on whether I and J are equal 
or not. This is the case if and only if the parameter matrices are of the form 

IJ = CT + {0 + -0~)5 IJ (15) 

with diagonal entries + and off-diagonal entries 0~ for O = P,Q, R, M, N. An explicit calculation 
shows that the matrix U simultaneously diagonalizes the parameter matrices, so that 

O = U-0-UT 1 = diag(Oi,Oo,...,Oo), (16) 
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where Oq = + — O and 0\ = + + (N — 1)0 are the eigenvalues of O. Further introducing 
t) = U_ ■ h ab) the most general linearized field equations can be written in the equivalent form 

= & ab = ^ ■ dPd (a*l b)p + Q ■ + 1 • d a dti + m ■ d p d\ q r, ab + S • airr, ab , (17) 

where the parameter matrices 21 are now diagonal matrices of the form (|16p . Thus, 

the equations decouple and we can write them as 

= = P x SPd {c fy\ )p + QiU^ ah + Rxdad^ 1 + M^Pffoab + NiO^Vab , (18a) 
= < 6 = P <9 p d (a rj* )p + Q n% 6 + i? £W + M Wr,^ + JVbDbSaft • (18b) 

A similar decomposition can be applied to equations Q, ((HI) and in terms of the quantities 
2 = U ■ I. These equations are invariant under gauge transformations if and only if they can be 
expressed in terms of gauge-invariant quantities only, i.e., they must not depend on the gauge- 
dependent quantities 3\, 3\ and 3'^. This is the case if and only if the eigenvalues of the parameter 
matrices satisfy the conditions 

Pi + 2Qi = Pi + 2Ry = Mi + Ni = , (19) 

as we have shown explicitly in [1] . In the following we will consider only multimetric theories which 
satisfy these gauge invariance conditions. 



B. Bianchi identities 



In the preceding subsection we have shown that the parameter matrices P,Q,R,^£,N_ are 
significantly restricted by the assumptions (0)-(^) stated in the introduction. In this subsection we 
will derive further restrictions which originate from the Bianchi identities. These are expected to 
hold since the theories we consider are derived from an action according to assumption @). Writing 
the gravitational part So[g , • • • , g ] of the action as an integral of the Lagrangian density C and 
demanding that it is invariant under diffeomorphisms generated by an arbitrary vector field £ leads 
to the condition 

o = 6tS Q = -2 1 d' x g 77vi \^^r\ ale ■ (20) 

Using the perturbation ansatz ((2|) and the basis transformation (|13[) this reduces to the linearized 
Bianchi identity 

= d a &\ h = (\p x + Q ± ) OdXb + (Vi + M t ) dbdP&lL. + (Pi + iVipdfef) 1 , (21) 
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which is a geometric identity. It then follows that the parameters must satisfy the additional 
constraints 

Pi = -2Qi = -2Mt , R ± = -Ni . (22) 

Together with the gauge invariance conditions ()19|) these constraints restrict the eigenvalues 
Pi,Qi,R%, Mi, Ni of the parameter matrices to only one free parameter, e.g., Pi, and the remain- 
ing parameters are determined as 

P 1 = -2Qi = -2Ri = -2Mi = 2Ni . (23) 

Note that in the single-metric case, in which the parameter matrices P_,Q,R,M_,N_ are replaced 
by their unique eigenvalues Pi, Qi, Ri, Mi, Ni, this determines the linearized field equations of a 
single-metric theory to be identical to those of general relativity, up to a constant factor. 

The situation is different for the remaining linear combinations f)* b of the metric perturbations. 
For these we do not obtain any constraints from the diffeomorphism-invariance of the gravitational 
action, since they are gauge invariants according to ()10p . However, we do obtain constraints from 
the fact that the matter action SmW ^ ,/ ] for each of the standard model copies is diffeomorphism 
invariant and so the corresponding energy- momentum tensors T^ b are conserved. Since the field 
equations are of the form K ab = 87rGNT ab according to assumption fzzj), this leads to the remaining 
linearized Bianchi identities 

= d a & ab = Qp + Q ) Vd% b + Qp + M ) d b d*>d% q + (R + NtiUdbtf , (24) 

which must be satisfied by all solutions of the linearized field equations ([3]). Note that in contrast 
to the Bianchi identity (|2ip , which follows from the diffeomorphism-invariance of the gravitational 
action, these additional Bianchi identities are in general not geometric identities, and are not 
necessarily satisfied by arbitrary metric perturbations. In order to satisfy (|24p we are thus left 
with two possibilities: 

(i) The Bianchi identities (|24p are geometric identities and satisfied by arbitrary perturbations 
t) l ab of the metric tensors. 

This is the case if and only if the eigenvalues of the parameter matrices satisfy the additional 
constraints 



P = -2Q = -2M , R Q = -N . 
It then follows that we are left with only three free parameters, e.g., Pi, Pq, Rq. 



(25) 
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(ii) The Bianchi identities (|24p are not geometric identities and satisfied only by solutions of the 
gravitational field equations. 

In this case we do not obtain any additional constraints on the parameter matrices. Note 
that even in this less restrictive case the Bianchi identities (|24p are implied by the grav- 
itational field equations, taking into account that the source of the gravitational field is 
given by the matter energy momentum tensors, which must be divergence-free due to the 
diffeomorphism-invariance of the matter action. Thus, (|24|) must be satisfied by all solutions 
of the gravitational field equations. 

In the remainder of this article we will not assume that the Bianchi conditions (|25p on the pa- 
rameters are satisfied in general, and instead use the smaller set (|23|) of combined Bianchi and 
gauge conditions. We will explicitly show how the different sets of conditions influence the possible 
wave-like solutions in the following section. 

C. Wave ansatz 

We are now in the position to explicitly construct wave-like solutions to the gauge invariant 
field equations ([7]) , (|SJ) and © . For simplicity, we apply the basis transformation (fl~3]) , so that all 
occurring parameter matrices become diagonal and the field equations are written in terms of the 
gauge-invariant quantities 3. For these quantities we use the wave ansatz 

3 = 3ex.p(ik a x a ) , (26) 

for a single Fourier mode, where 3 are constants which determine the amplitude of the wave and 
k a is a constant covector. The aim of this section is to compute for which amplitudes and wave 
covectors the field equations are satisfied. 

First we solve the tensor equations ([9]) . Using the basis transformation (|13p , these take the form 

= 2Q l D3 1 a p , = 2Q n3 i a(5 . (27) 

If one of the eigenvalues Qo,Qi vanishes, the corresponding equation is satisfied identically. In 
this case the linearized field equations ((3|) are not sufficient to solve for some of the modes 3^ and 
a calculation based on the full nonlinear field equations is required. We will not attempt such a 
calculation in this article and restrict ourselves to the case that both Qq and Q\ are nonzero. We 
then insert the wave ansatz (|26p and obtain the equations 

= -2Qik a k a 3^ , = -2Q k a k a 3^ . (28) 
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We immediately see that non- vanishing wave-like solutions exist if and only if the wave covector 
k a is null, k a k a = 0, while the amplitudes 3^ may be arbitrary. 

Second we consider the vector equations (JSj) . After applying the basis transformation and 
inserting the wave ansatz, these can be written in the form 







y 1 



-fPi^o 
\P kl-Q Q k a k a -i(p + 2Q )k k a k a 
-±(P + 2Q )k a k a 



iPo&Q 



(29a) 
(29b) 



where we already used the gauge invariance conditions (|19p to eliminate Q±. Consequently, the 
pure gauge quantity 3'^ drops out and the field equations depend only on the physical degrees of 
freedom. From equation (J29aj we see that we can solve for the quantities 3"^ using the linearized 
field equations only if Pi ^ 0. It then follows that there are no wave-like solutions for 3**, since 
the field equations require that both fco an d k a must vanish. Similarly, a solvable equation for 3 l a 
and 3'a requires that both Q$ and Pq + 2Qq are nonzero. It then follows that wave-like solutions 
exist if and only if the determinant of the matrix in equation (|29bj) vanishes, 



-MP + 2Q )(k a k 



a\2 



o, 



(30) 



i.e., for null waves, k a k a = 0. For these solutions 3'^ is always allowed to be nonzero, while 3 l a 
is allowed to be nonzero only if Pq = 0. Note that Po + 2Qq for theories in which the Bianchi 
identities (|24p are geometric identities and in which the parameters satisfy the conditions ()25|) . In 
this case the linearized field equations are not solvable for the vector potential 3^. The remaining 
vector potential 3 l a must vanish for Pq ^ and cannot be determined from the linearized equations 
for P = 0. 

Finally we discuss the scalar equations ©. We can proceed in complete analogy to the vector 
equations shown above. We apply the basis transformation (|13p . insert the wave ansatz (|26p and 
make use of the gauge invariance conditions (|19p in order to eliminate the parameters Qi, R\ and 
N\. The equations for the first component J 1 can then be written in matrix form as 



/ -(P x + 2M 1 )k a k a AM x k a k a - 3(Pi + 2M l )kl o\ 
2iP 1 k 

(Pi + QM{)k a k a -4(Pi + 3Mi)k a k a + 3(P X + 6Mi)^ 
-Pi Pi 0/ 



~4 



(31) 
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We immediately see that these do not depend on the pure gauge quantities 3% and 3\ as a conse- 
quence of the gauge invariance conditions. For the gauge-invariant quantities 3\ and 3^ we m ust 
distinguish two different cases. If Pi 7^ 0, it follows from the second component equation of (|3ip 
that there are no wave-like solutions for 3\. From the last component equation we further see that 
3\ = 3 2 and thus there are no wave- like solutions for 3\ either. In the case Pi = the equations 
are not sufficient to determine the quantities 3\ and 3\. 

For the remaining quantities 3 % we proceed similarly and write the field equations in matrix 
form. As we already have done for the vector equation (|29bp . we calculate the determinant of the 
occurring matrix, which takes the form 

{3Qo(Po + 2Q ) [Qo(Po + Qo + Ro) + M (Q - 3P ) + N (3P + 4Q )] }(k a k a f = . (32) 

Again we distinguish two cases. If the constant factor in curly brackets vanishes, the linearized 
field equations are not sufficient to solve for the quantities 3 l . Otherwise, the field equations can 
be solved by the wave ansatz if and only if the determinant vanishes, which is the case for null 
waves. A quick calculation shows that the solutions take the form 



Note further that if the parameters satisfy the Bianchi conditions (|25p . the determinant (|32p 
vanishes identically so that the linearized field equations are not solvable for the scalar potentials 3 l . 

This result completes our discussion of gravitational waves in the gauge-invariant formalism. 
We have shown that we need to impose several conditions on the eigenvalues of the parameter 
matrices P, Q,R, M_, N_. From assumption Ud\l we concluded that they must be of the form (|15p 
and can be diagonalized according to equation (fTU|) . The gauge conditions (fTUj) guarantee that 
the linearized field equations are gauge-invariant and thus depend only on the physical degrees of 
freedom. The solvability conditions 

Qo^Q, Pi = -2Qi/0, P + 2Qo/0, 
Qo(Po + Qo + Ro) + M (Qo - 3R ) + N (3P + 4Q ) + (34) 

allow a complete treatment of gravitational waves using the linearized field equations. Under these 
conditions, wave-like solutions for the scalar, vector and tensor components of the gauge-invariant 
quantities exist if and only if the wave covector k a is null, k a k a = 0. Since we treat gravitational 
waves as a small perturbation of the flat Minkowski background, whose null directions govern the 
propagation of light within our approximation, it then follows that the speed of gravitational waves 
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equals the speed of light, v g = c. In the more general case that a concrete theory does not satisfy 
all of the solvability conditions (|34p . a calculation based on the full non-linear field equations is 
necessary to determine whether additional wave-like solutions exist. Since we do not perform such 
a calculation in this article, we restrict ourselves to the null wave solutions we have found so far 
and determine their possible polarizations in the following section. 



III. NEWMAN-PENROSE FORMALISM AND POSSIBLE POLARIZATIONS 



In the preceding section we have shown that the most general linearized vacuum field equa- 
tions ([3]) can be solved by the wave ansatz (|26p only if the wave covector is null. We will now 
turn our focus to the possible polarizations of gravitational waves. Since we are dealing only 



enrose 
in 



with null waves, the polarizations can easily be decomposed by employing the Newman- 
formalism introduced in [23]. We will then employ the classification scheme detailed in 
order to determine the E(2) class of multimetric gravity, which could be measured by the upcom- 
ing gravitational wave experiments. The connection between these experiments and the class of 
multimetric theories we consider in this article is established by assumption @) stated in the intro- 
duction. It follows from this assumption that a gravitational wave experiment built up from visible 
matter, i.e., only one standard model copy ^f 1 , is sensitive only to the corresponding metric g\ b . 
Thus, it is sufficient to determine the possible polarizations of wave-like solutions for the metric 
g\ h . From assumption fwj) it follows further that the possible polarizations are the same for all 
sectors (^ / ,<7 / ). 

Basic ingredient of the Newman-Penrose formalism is a convenient double null basis of the 
tangent space. In the following, we will use the notation of [14] and denote the basis vectors by 
l a , n a , m a , rh a . In the x°,x a basis they take the form 



l a = (1,0, 0,1) 



1 



n 



(1,0,0,-1) 



1 



m 



(o,m,o) 



l 



m 



V2 



(0,l,-i,0). (35) 



2 v-,-,-, - ^ 

In the new basis the flat Minkowski metric takes the form 

^ -10 0^ 
-10 

Vab = 

1 
\ 10/ 

We now consider a plane wave propagating in the positive x s direction. The wave covector then 



(36) 
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takes the form k a = —ujl a and the metric perturbations are given by 

ha b = ha»f* P -* ) = h b j" U (37) 

for the retarded time u = x° — x 3 . The effect of this wave on a set of test masses consisted by 
one type ^t 1 of standard model matter depends only on the Riemann tensor of the corresponding 
metric g\ h . As shown in |9(, the Riemann tensor of a plane wave is determined completely by the 
six so-called electric components. For the wave (|37p . these can be written as 



^.2 — Q^nini ~ ^2~^ ' ~^~ 3 — 2— nlnm ~ 2— n!,n ' m ~ 4— lm ~ ' 

lv lv _ _ 1 



^4 Mhnmnm I^nmnm ^— Trim ^— mm 1 —22 R-nmnm ^—mm ' (3^) 



where dots denote derivatives with respect to u. We now examine which of the components (|38j) 
may occur for gravitational waves satisfying the linearized field equations ([3]). Inserting the wave 
ansatz (|37|) we immediately see that the terms containing Q and N_ drop out, since Oh ab = for 
a null wave. Writing the curvature tensor K_ ab in the Newman-Penrose basis, we find that the five 
component equations 

= Kjl = K.mm = M-mfh = ¥±im = Kjfh (39) 

are satisfied identically, while the remaining five component equations take the form 

= K nn = 2E-h mm -(E + 2R)-h ln , (40a) 
= K in = ~(P + 2M)-h Ui (40b) 

= Knm = -^E-hlm, ( 40c ) 

= K nin = -^E-h lrh , (40d) 
= K m m=K-ki- (40e) 

Recall that the parameter matrices can be brought into diagonal form using the basis transforma- 
tion (fl~6j) so that the field equations decouple as shown in (fTBj) . Applying this decomposition to 
equation (|40ap we obtain 

= & nn = 2R X ■ blrn ~ (Pi + Ml) ■ bi„ , (41a) 
= 4 = 2R ■ ~if mih - (P + 2Rq) • b?„ , (41b) 

and similarly for the remaining four component equations. Let h a b denote one of the metric 
perturbations fj^ and P, R, M the corresponding eigenvalues of the parameter matrices P. R, M_. 
We distinguish the following cases: 
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• M = P = 0: In this case equations (|40bp and (|40ep are satisfied identically for arbitrary 
amplitudes hu. For waves of this type the corresponding component R n ini = — 6\l/oof the 
Riemann tensor is allowed to be nonzero. Following the classification detailed in [9|, they 
belong to the E(2) class IL3. 

• M / and P = 0: Equation (|40ep forbids waves with a nonzero amplitude hu, and thus 
^2 = 0. Equations ()40cp and (|40dp are satisfied identically for arbitrary amplitudes h\ m and 
^/flj. It then follows that R n inm = is allowed to be nonzero. Waves of this type belong 
to the E(2) class III 5 . 

• P ^ and P + 2R / 0: For P / it follows from equations (t40b|) . ([TO . (|40dl and (|4HeD 
that hu = hi m = him = 0, and thus ^2 = ^3 = 0. The remaining equation (|40ap is solved 
for 2Rh m fh = (P + 2R)h>in- Hence, the corresponding component R n mnm = — ^22 of the 
Riemann tensor is allowed to be nonzero and the wave belongs to E(2) class N3. 

• P = —2R ^ 0: This is the most restrictive case. Equations (|4Up are satisfied only for 
hu = h[ m = hi m = h m rh = 0, and thus ^2 = ^3 = ^22 = 0. The only allowed polarization is 
Rnmnm = — ^4 and the wave belongs to E(2) class N2. 

The classification can be summarized in a convenient graphical form. The following diagram shows 
how the E(2) class is determined by the values of the parameters P, R, M in the linearized field 
equations: 

P 



P + 2R 

N 2 N 3 III5 H 6 

2 tensors +1 scalar +2 vectors +1 scalar 

The parameters P,R,M in this diagram are either the set P\, R\, M\, which yields the E(2) class 
for waves of type \) l ab , or Pq, Rq, Mq, which instead yields the E(2) class for waves of type t) l ab . We 
thus obtain two E(2) classes for the different linear combinations of the metric perturbations h^. 

We finally turn our focus to the viewpoint of a physical observer. From assumption fj) in the 
introduction it follows that any experimental setup consisting of visible matter only, i.e., of only 
one copy ^ x of the standard model, is affected by only one metric tensor g x ah . A gravitational 
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wave detector consisting of visible matter can thus measure only one of the Riemann tensors 
R^ bcd - From the basis transformation (|13p it follows further that each Riemann tensor R^ bcd 
depends on all metric perturbations \) T ab . As a consequence, it is not possible to measure the metric 
perturbations rj^ 6 separately. A gravitational wave experiment can only indicate whether a wave 
with certain polarization exists for any of the linear combinations fr 6 . Thus, only the larger of the 
two aforementioned E(2) classes can be determined. 



IV. EXAMPLES 

In the previous section UlTl we constructed a formalism to calculate the E(2) class of multimetric 
gravity theories that determines the possible polarizations of gravitational waves. It turned out 
that the E(2) class is fully determined by the eigenvalues of the parameter matrices P, Q, R, M, N 
in the linearized vacuum field equations ([3]). We will now apply this classification to a number of 
example theories. 



A. General relativity 

Although general relativity is not a multimetric theory, we can apply a slightly modified version 
of the calculations presented in this article. For the case of N = 1 metric tensors, we replace 
the parameter matrices in the linearized vacuum field equations by their unique eigenvalues. For 
general relativity, these take the values 

Pl = l, Q 1 =R 1 =M 1 = ~, JVi = i. (42) 

One easily checks that they satisfy the gauge invariance conditions (|19p and the Bianchi condi- 
tions (|22p . Following the calculation presented in subsection III CI one finds that the wave solutions 
are completely determined by the linearized field equations, and that the only permitted non-zero 
amplitude is the tensor jLg. Finally, a comparison of the parameters (|42p with the diagram at the 
end of section [III] correctly shows that the E(2) class of general relativity is N2. Note that this is 
the generic case for single-metric theories due to the parameter constraints (|23|) . 
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B. A simple multimetric theory 

A simple class of multimetric gravity theories with N > 2 metrics, which also contains the 
theories presented in 0, [^j] , is given by the gravitational action 

1 r r N N ( 

sob 1 , ...,g N ] = — d^x^-o x !i" J "'u + Y.a IlJ yrfu 

n J I I,J=1 1=1 \ 

i el i „Z'l c' ,.,ci A- cl "' i i<i .J I ,. el i» c' " i ,, Z'l in c'ln \ \ (43) 



+ uSUS'j + vS'kS 1 % + wS 1 fc im 5 J m jk + g 1 Vmn rS 1 m lk S J n jl + sS 1 m ij S 1 n kl 



where the connection difference tensors S Ik ij, S \ are defined as 

1 N 

qIJ i -pi i -pJ i qIJ qIJ k oJ i \ oUi qJ a J k (AA\ 

D jk — L jk 1 jk j <-> j ~ D jk > »-> jk — / 4 D jk j D j ~ D jk > v") 

1=1 

the volume form is given by go = II/Li {d 1 ) 1 ^ '> an d x,y,u,v,w,r,s are constant parameters. 
Starting from the action (|43p . we derive the gravitational field equations by variation with respect 
to the metric tensors and use the perturbation ansatz ([2]) to keep only the terms of linear order 0(h). 
This yields the eigenvalues of the parameter matrices 

Nx - v + 2s 

Pi = -2Qi = -2R X = -2M X = 2Nx = Nx + y, R = M = , 

^ Nx — y + w — 3r , T — Nx — y — u + v — s 
P = - Nx + y - w + r-2s, Q = U — , N = y — , (45) 

which satisfy the gauge invariance conditions (|19p and the Bianchi conditions ()22p since we started 
from a diffeomorphism-invariant action. The extended Bianchi conditions (|25|) are satisfied if and 
only if the parameters satisfy the constraints 

Q=s+r=u+v+w=u+y— s . (46) 

It thus follows from our discussion of the Bianchi identities in subsection III Bl that in the generic 
case, in which the conditions (|46p are not satisfied, the Bianchi identities (|24p are not geometric 
identities, but satisfied only by solutions of the gravitational field equations. Next, we apply the 
linearized multimetric extension of the parametrized post-Newtonian formalism detailed in Q]. 
Consistency with solar system measurements of the PPN parameters requires 

1 6-JV Ar 6-N Ar 1 

y = j—j^-Nx, v = — —-Nx + 2u, w = - 4 _ 2Ar + Nx-3u, r = -j—^ + Nx-u, 

(47) 
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which leaves us with only three free parameters x, u, s and restricts the number of metrics to N > 2. 
In terms of these remaining parameters the eigenvalues of the parameter matrices take the values 

Pi = -2Qi = -2R 1 = -2Mi = 2iVi = , R = M = - £~ m +Nx-u + s, 

„ 6-iV AT „ 14-iV -Nx + u-s /A , 

Po = I ^-2Nx + 2u-2s, Q = --, N = ^-^ + . (48) 

Using the calculation of subsection III CI we find the following wave solutions for the gauge invariant 
potentials 3: 

• Tensor modes: 

Since both Q\ and Qo are nonzero, wave- like solutions for all tensor potentials 3^ exist. 

• Vector modes: 

Prom Pi 7^ it follows that there are no wave-like solutions for the vector potential 3*. If 
the parameters satisfy the condition 

Nx-u + s = , (49) 

the linearized field equations are not sufficient to determine the vector potentials 3„ an d the 
remaining vector potentials 3 l a must vanish. Otherwise, wave-like solutions for the vector 
potentials 3'^ exist, and wave-like solutions for 3 l a exist if and only if the parameters satisfy 

6 — N 

Nx-u + s = . (50) 

8 - AN v ' 

• Scalar modes: 

Prom Pi ^ it follows that there are no wave- like solutions for the scalar potentials J\, 3\. If 
the parameters satisfy the conditions (|49p , the linearized field equations are not sufficient to 
determine the scalar potentials 3\, 3 l 2 , J|, 3\. Otherwise, wave-like solutions for 3\,3 l 2 ,3 l 3 ,3\ 
exist and satisfy (|33|) . 

Finally, we determine the E(2) class of our example theory. A comparison of the eigenvalues (|48p 
with the classification detailed in section IIIII shows that the E(2) class for the symmetric linear 
combinations f)^ 6 of metric perturbations, for which the eigenvalues Pi, Pi, Mi of the parameter 
matrices are relevant, is N2. The E(2) class for the remaining linear combinations t) l b , and thus 
the effective E(2) class of the theory, is IL3 if the parameters satisfy the condition (j50j) . and N2 
otherwise. This means that the only polarizations that can be measured by a gravitational wave 
experiment are the two tensor polarizations which are also present in general relativity, unless the 
parameters satisfy ()5Q|) . in which case all six possible polarizations of gravitational waves may be 
present. 
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V. CONCLUSION 

In this article we have discussed the propagation of gravitational waves in theories with N > 2 
metric tensors and a corresponding number of standard model copies. These theories were designed 
to explain the observed cosmological late-time acceleration while being consistent with solar system 
experiments. We have examined two characteristic properties of gravitational waves: their propa- 
gation velocity and their polarization. We have shown that in a weak-field approximation around 
a flat, maximally symmetric Minkowski background all gravitational waves propagate at the speed 
of light. Using the Newman-Penrose formalism we found that there are always two tensor modes; 
in addition two vector modes and two scalar modes may exist. In terms of E(2) representations 
this means that multimetric gravity theories can be of class N2, N3, III5 or II^. 

We then applied our construction to two examples. First, we discussed general relativity and 
showed that our formalism can also be applied to the special case N = 1 of a single-metric gravity 
theory. We re-obtained the well-known result that general relativity is of class N2, i.e., there are 
only two tensor polarizations of gravitational waves. Second, we applied our construction to a class 
of multimetric gravity theories including the theories proposed in jj, [J and showed that these are 
either of class N2 or IL3, depending on the choice of parameters. 

Our results connect the theoretical framework of multimetric gravity to the physics of gravita- 
tional waves, which is the subject of several current and upcoming experiments. It is expected that 
these will be able to measure both the propagation velocity and the polarization of gravitational 
waves, which are the two properties addressed in this article. A question of particular interest 
is whether scalar or vector polarizations will be detected. While these do not exist in general 
relativity, they are allowed in certain multimetric gravity theories. 

Now that we examined the propagation of gravitational waves in multimetric gravity theories, 
the next task that should be performed is to discuss their production by sources such as binary 
stars. Further research on this topic should show which of the propagating wave polarizations are 
emitted from a given source. Moreover, quantitative calculations should yield the amplitude of the 
emitted waves in the different metric sectors. Since energy may be emitted in all metric sectors, but 
only one of them is visible to gravitational wave detectors, one might expect a difference between 
the directly observed energy emission and the total energy loss inferred from the orbital decay. 
This could provide another testbed for multimetric gravity both by the upcoming gravitational 



wave experiments 



lfl 



ll| | and existing observations of binary pulsars 2514271] . 



Finally, it should be examined how the presence of cosmic gravitational fields affects both 
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electromagnetic and gravitational radiation as they propagate from a common source, such as a 
supernova or a binary pulsar, towards our solar system. Calculations of this type are crucial for 
the interpretation of experiments which compare the arrival times of both types of radiation, as 
they might undergo a different Shapiro delay [if], [n| or a different gravitational lensing. 
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